Competition of crystal field splitting and Hund's rule coupling in two-orbital magnetic 

metal-insulator transitions 
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Competition of crystal field splitting and Hund's rule coupling in magnetic metal-insulator transitions of half- 
filled two-orbital Hubbard model is investigated by multi-orbital slave-boson mean field theory. We show that 
with the increase of Coulomb correlation, the system firstly transits from a paramagnetic (PM) metal to a Neel 
antiferromagnetic (AFM) Mott insulator, or a nonmagnetic orbital insulator, depending on the competition of 
crystal field splitting and the Hund's rule coupling. The different AFM Mott insulator, PM metal and orbital 
insulating phase are none, partially and fully orbital polarized, respectively. For a small J H and a finite crystal 
field, the orbital insulator is robust. Although the system is nonmagnetic, the phase boundary of the orbital 
insulator transition obviously shifts to the small U regime after the magnetic correlations is taken into account. 
These results demonstrate that large crystal field splitting favors the formation of the orbital insulating phase, 
while large Hund's rule coupling tends to destroy it, driving the low-spin to high-spin transition. 

PACS numbers: 71.30. +h, 75.30.Kz, 71.10.Hf, 71.27.+a, 71.10.Fd 



I. INTRODUCTION 

Mott-Hubbard metal-insulator transitions (MIT) in multi- 
orbital Hubbard models have been extensively studied since 
there are rich phase diagrams '~ 12 . A few of many-body ap- 
proaches have been employed to study the paramagnetic Mott 
transition, such as dynamic mean field theory (DMFT) 1 -^ and 
Kotliar-Ruckenstein slave-boson method 910 , etc. With the 
slave spin technique, Yu and Si recently showed that a para- 
magnetic orbital insulating phase can exist when the Hund's 
coupling approaches zero 11 . However, generally speaking, 
MIT transitions are usually accompanied by magnetic tran- 
sitions. How spin correlation affects the nature of MIT is sel- 
dom discussed. 

On the other hand, the crystal field splitting plays an im- 
portant role in the properties of multi-orbital systems. It is 
not clear whether the interplay between crystal field splitting 
and spin exchange splitting is a key factor for antiferromag- 
netic (AFM) MIT. Though Hasegawa 13 studied the spin cor- 
relation effect in half-filled doubly degenerate Hubbard model 
and found the Hund's coupling plays a very considerable role 
in magnetization, he did not find the MIT at half filling. It 
deserves to explore the influence of crystal field splitting on 
the MIT in the presence of magnetism. And how the crystal 
field and spin correlation modify the orbital insulating phase 
is also an interesting issue. 

In this paper, the multi-orbital Kotliar-Ruckenstein slave- 
boson method is generalized with spin degree of freedom. 
The effects of crystal field splitting and Hund's rule coupling 
on Mott transition are also discussed. We find that besides 



Jh = 0, orbital insulator can exist for small enough but fi- 
nite Jh, different from Yu and Si's results 11 . The magnetic 
phase diagram as a function of U and crystal field splitting 
shows that the phase boundary can be seriously changed by 
the Hund's rule coupling and crystal field splitting, demon- 
strating the contrary role of these two factors. The rest of this 
paper is organized as follows: in Sec. II we describe the model 
Hamiltonian and theoretical approach. The effects of crystal 
field splitting are discussed and corresponding phase diagram 
are presented in Sec. III. The final section is devoted to the 
conclusion remarks. 



II. MODEL HAMILTONIAN AND METHODS 

The two-orbital Hubbard model adopted in this paper is as 
follows: 



H = H + Hi 



(1) 
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where c f is a creation operator of an electron with the or- 
bital index a and spin <x at the lattice site i , tii aa - is the cor- 
responding occupation number operator. The hopping inte- 
gral between orbitals a and /3 is denoted by t a p. The intra- 
band (inter-band) Coulomb repulsion and Hund's rule cou- 
pling are denoted by U (U' ) and Jh, respectively. Here 
we set U' — U — 2Jh- We introduce new boson operators 
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e, p, d, b, t and q, which denote the empty, single occupation, 
double occupation in two orbital, double occupation in one 
orbital, triple occupation, and quadruple occupation, respec- 
tively. The physical electron creation operator cl is repre- 
sented as Z i(Hr /. , where f! is quasi-particle creation oper- 
ator. The renormalization factor Z iaiT is given by 



= Q 2 



+ 2 icr-^^ + 



(l - G,^)- 



(4) 



where 



Qiao- - p\ a(T Piaa- + b\ a b ia + ^ d^.d^^ + ^ 



(5) 



In the present slave boson states, the normalization con- 
straint and Fermi number constraint maintaining physical 
Hilbert space are as follows: 



1 = e\e 
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Our numerical method which is used to search the minima 
energy is based on the pattern search method, the gradient 
method and the Rosenbrock method. The normalization con- 
straint given by Eq.© must be held at all time in our numer- 
ical calculation. The Fermion number constraint is enforced 
by using the penalty function method. Our numerical calcula- 
tions are performed for simple square lattice with the nearest- 
neighbor hopping. Throughout this paper we set the ratio of 
two intraorbital hopping integrals as ?22/fn = 0.5, and the av- 
erage electron number per site is 2. We measure the energies 
in units of the bandwidth of the orbit- 1, 2D \ = 8?n. 



HI. NUMERICAL RESULTS 

We first present the magnetic phase diagram of two-orbital 
Hubbard model at half filling as the functions of U and A in 
Fig. 1. For a large Hund's coupling, e.g. Jh = 0.2(7, there are 
only paramagnetic (PM) metallic phase when U < U c \ and 
AFM insulating phase for U > U c \ in the U - A phase dia- 
gram. We notice that the magnetic MIT phase boundary U c \ 
shifts downward with the increase of the crystal field splitting. 
When the system transits from a PM metal to an AFM insu- 
lator, a fraction of electrons transfer from a lower occupied 
orbit to an upper unoccupied orbit. Consequently, the spin 
exchange splitting should overcome the crystal field splitting, 
leading to the decrease in the groundstate energy. We expect 
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FIG. 1: Magnetic phase diagram in the U - A plane with tn = 0, 
where U c i and U C 2 are the PM metal - AFM insulator boundary at 
Jh = 0.2U and J H = 0.02(7 respectively. U c j and U C 4 are the orbital 
insulating phase boundaries at J H = 0.02(7 with and without spin 
degree of freedom respectively. 



that the magnetic MIT is first order, since it occurs only when 
the spin exchange splitting is superior to the crystal field split- 
ting. The orbital polarization induced by the crystal field split- 
ting is also expected to decrease sharply in the AFM insulating 
phase. 

For the small Hund's rule coupling situation with Jh = 

0. 02 U, the magnetic phase diagram is richer than that with 
J H = 0.2U. Besides the PM metallic and AFM insulating 
phases, an orbital insulating state is stable when U > U c i, 

1. e. the system becomes completely orbital polarized with- 
out magnetism in the large A region. As shown in FigQ] the 
present PM-AFM phase boundary U& shifts to large U, in 
comparison with the phase boundary U c \ with Jh = Q.2U. In 
the comparable large U and A and small Jh region, two elec- 
trons fall into the lower energy orbit. As a result, the system 
undergoes an AFM/PM to orbital insulator transition. 

The orbital insulating phase boundary without spin corre- 
lation is also plotted as U C 4 in Fig[T] which is in agreement 
with Werner and Millis' result 4 . Obviously, the present criti- 
cal value for orbital insulator, U C 3, is considerably lower than 
that without spin correlation, U C 4, though there is no sublattice 
magnetic moment in the orbital insulating phase. The reason 
is that the modulation of spin correlation on electronic spec- 
trum is in favor of the occurrence of the orbital insulator. 

To understand the electron distributions more clearly, we 
present the dependence of the boson occupancy probabilities 
as the functions of U for different crystal field splittings in 
Figf2] In the case of Jh = 0.02 U, the spin-singlet double 
occupation probability Z?2 is dominant over the whole corre- 
lation range when the crystal field splitting A is comparably 
large. In the case of Jh = 0.217, the spin-triplet double occu- 
pation probability d mT is dominant, as seen in FigfJ] Figf2ta) 
shows that when the orbital splitting A = 0, the spin-up and 
spin-down occupancies split and the ground state of the sys- 
tem is AFM at finite U. From the inset of Fig|2ja), one no- 
tices that dn considerably rises with the increase of U, since 
the large U and Jh are in favor of the electron distribution on 
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FIG. 2: Boson occupancy probabilities of sublattice as functions of 
U for different crystal field splitting at n = 2. Theoretical parameters: 
J H = 0.2(/, tu = 0.2f u . (a) A=0; (b) A=0.1. Insets in (a) and (b) 
show the full scale of the dependence of occupancies on U 



different orbits with parallel spin, resulting in the dominantly 
increase of d^. And the evolution of various orbital occupan- 
cies probabilities on electron correlation qualitatively agrees 
with Hasegawa's results^. 

However, in the presence of crystal field splitting, the bo- 
son occupancy probabilities shown in Figj2f b) are very differ- 
ent from Fig|2a). With increasing U, the system undergoes a 
first-order MIT transition and transits from PM to AFM phase 
at U — U c \ since spin exchange splitting becomes larger than 
crystal field splitting and spin gap comes into being. At the 
same time, the spin-up and spin-down occupancies do not split 
until U > U c \. When U > U c \, various occupancies mono- 
tonically decrease with increasing U, except for dn, which 
steeply increases and gradually approaches a saturation, as 
one can see in the inset of Fig|2jb). 

The evolution of the particle distributions in two orbits is 
plotted in Fig[3] Figj3ja) shows that the sublattice magnetic 
moment contributed from narrow band is lager than that from 
the wide band, since the spin exchange splitting is stronger in 
the narrow band. In the present two-orbital system with A = 
and large Jh in Figf3ja), the orbital polarization is always sup- 
pressed since the gravity centers of two bands do not shift 
with each other, the orbital singlet occupation dominates both 
in the PM metallic and in the AFM insulating phases. And 
from Figj3j), it is found that although the crystal field split- 
ting induces the partial orbital polarization in the PM metallic 
phase, it leads to completely unpolarized in the AFM insulat- 
ing phase since the each localized electron is allocated each 
orbits, eliminating the orbital polarization. 

The correlation dependence of the particle distributions and 
the renormalizations factors of two orbits at Jh = 0.02(7 and 
A = 0.18 are plotted in FigUa)and (c). The corresponding 
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FIG. 3: Evolution of Electron number distributions in two orbits, (a) 
A = and (b) A = 0. 1. Other parameters are the same to Figj2] 
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FIG. 4: Dependence of electron occupation numbers and band renor- 
malization factors of two orbits on Coulomb correlation, (a) and (c) 
with spin correlation, (b) and (d) without spin correlation. Theoreti- 
cal parameters: A = 0.18, J H = 0.02(7, t l2 = 0. 



results without considering spin exchange splitting are also 
plotted in FigHtb)and (d) for comparison. Both Figj4ja) and 
Hfb) show that the difference of particle numbers in two or- 
bits, i.e. orbital polarization, increases with the increase of U, 
and the system becomes completely orbital polarized in the 
orbital insulating phase when U > U c j, accompanied with the 
PM-AFM MIT. The corresponding renormalization factors in 
FigHfc) and|4|d) show that the MIT simultaneously occur in 
the two orbits, even if in the presence of finite crystal field 
splitting. Comparing FigHfa) with Figj4fb) and FigHfc) with 
FigHfd), we find that the spin correlation greatly reduces the 
critical value U c of the MIT, therefore favors the occurrence 
of the MIT. 

Our numerical results for different inter-orbital hopping 
do not qualitatively affect the magnetic MIT, only the phase 
boundaries show a slight shift. We find that different from the 
statement by Yu and Si 11 , in addition to Jh = 0, the orbital 
insulating phase can exist for finite but small Jh, since the 
competition between crystal field splitting and the Hund's rule 
coupling may stabilize the orbital insulating phase. It is in- 
teresting to expect that considerable lattice distortion may be 
associated with the occurrence of orbital insulator due to the 
first-order phase transition with large variation in orbital po- 
larization. Further, as seen in FigJT] |2]and|3] the high-spin to 
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low-spin state transition is clearly found in our present study, 
proving Millis et a/.'s conjecture in their paramagnetic DMFT 
study 4 , which indicates that the MIT with magnetism is much 
more rich than that without magnetism. 

IV. SUMMARY AND CONCLUSION 

We have performed a comparative study on the alterna- 
tive effects of crystal field splitting and Hund's rule cou- 
pling on MIT of half-filled two-orbital asymmetric Hubbard 
model with and without magnetism by means of the Kotliar- 
Ruckenstein slave boson approach. Our results show that the 
MIT is first order accompanied with a magnetic transition in 
the presence of crystal field splitting. The critical value U c for 
MIT considerably reduces after taking into account the spin 
correlation. Meanwhile, the orbital insulator can exist in the 



system with large crystal field splitting and small Hund's rule 
coupling. With the increase of crystal field splitting, the sys- 
tem may enter from a high-spin AFM state to a low-spin or- 
bital insulating phase in the small Jh region. These results 
demonstrate that the competition between the crystal field 
splitting and the Hund's rule coupling Jh is very important 
for the metal-insulator transition with spin degree of freedom. 
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